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Abstract 

We discuss an algebraic treatment of three-body systems in terms of a U{7) spectrum generating 
algebra. In particular, we develop the formalism for nonlinear configurations and present an algebraic 
description of vibrational and rotational excitations of symmetric (X3) and asymmetric tops (XY2 
and XYZ). The relevant point group symmetry is incorporated exactly. 
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1 Introduction 



The study of few-body problems has played an important role in many fields of physics and chemistry 
0. Over the years accurate methods have been developed to solve the few-body equations. The degree 
of sophistication required depends on the physical system, i.e. to solve the few-body problem in atomic 
physics requires a far higher accuracy than in hadronic physics. 

In recent years, the development and application of algebraic methods to the many-body problem {e.g. 
collective excitations in nuclei and molecules [||) has received considerable attention. In spectroscopic 
studies these algebraic methods provide a powerful tool to study symmetries and selection rules, to 
classify the basis states, and to calculate matrix elements. In this paper we discuss the application of 
algebraic methods to the few-body problem. Especially in the area of hadronic physics, which is that of 
strong interactions at low energies, for which exact solutions of QCD are unavailable, these methods may 
become very useful. The approach is based on the general criterion Q] to consider U{k + 1) as a spectrum 
generating algebra for a bound-state problem with k degrees of freedom, and assigning all states to the 
symmetric representation [N] of U{k + 1). For the fc = 5 quadrupole degrees of freedom in collective 
nuclei this led to the introduction of the U (6) interacting boson model j|] and for the A: = 3 dipole degrees 
of freedom in diatomic molecules to the U (4) vibron model § . 

In a three-body system the dynamics is determined by the six degrees of freedom of two relative 
vectors, which in the algebraic approach leads to a U (7) spectrum generating algebra. This model was 
originally introduced for a system of three identical objects ^, In this paper, we develop the general 
formalism for a system of three objects and discuss, in particular, an algebraic description of vibrational 
and rotational excitations of symmetric and asymmetric tops. Applications can be found in molecular 
physics (XYZ, XY2 and X3 molecules) and hadronic physics {qqq baryons). 



2 Algebraic treatment of a three-body system 

The internal motion of a three-body system can be described in terms of the relative Jacobi coordinates, 
p and A, which we choose as 

A = , \ = [wifi -I- 7712^*2 - {mi + 7712)^5] • (1) 

^ym^+ 1712 + ("^1 + ™2) 

Here rrii and r^; denote the mass and coordinate of the i-th object. For three identical objects with 
equal masses, Eq. (|l|) reduces to the Jacobi coordinates used in Instead of a formulation in terms 
of coordinates and momenta, we use the method of bosonic quantization in which we introduce a dipole 
boson with angular momentum and parity = 1^ for each independent relative coordinate, and an 
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auxiliary scalar boson with = 0+ 

pL-pU'^^ (m = -1,0,1). (2) 

The scalar boson does not represent an independent degree of freedom, but is added under the restriction 
that the total number of bosons N — Up -\- n\ + Ug is conserved. This procedure leads to a spectrum 
generating algebra of U{7) |^. For a system of interacting bosons the model space is spanned by 
the symmetric irreducible representation [N] of U{7), which contains the harmonic oscillator shells with 
n = Up + n\ — 0,1, . . . , N. The value of N determines the size of the model space. 

In case two or all three objects are identical, the Hamiltonian has to be invariant under the permutation 
group, iS'2 or 5*3, respectively. In the former case, in which we label the identical objects by 1 and 2, the 
permutation symmetry is determined by the transposition P(12), and in the latter case by both -P(12) 
and the cyclic permutation P(123) Q]. All other permutations can be expressed in terms of these two 
elementary ones. Algebraically, these operators can be transcribed as 



P(12) — expl—infip} 
27r 

T 



P(123) = cxp I ~i^F2} , (3) 



with 



F2 = -i 



J2 {plrnP\^n - pi^^Pp,,n) ■ (4) 



For three identical objects the eigenstates of the Hamiltonian are characterized by the irreducible repre- 
sentations of the S3 permutation group |0, However, in anticipation of the geometric analysis of the 
next sections, we use a labeling under the point group (which is isomorphic to S3): namely Ai and 
A2 for the one-dimensional symmetric and antisymmetric representations, and E for the two-dimensional 
mixed symmetry representation. The scalar boson of Eq. transforms as the symmetric represen- 
tation Ai , whereas the two dipole bosons pj^ ^ and p| ^ transform as the two components of the mixed 
symmetry representation, Ep and E\, respectively. For two identical objects the Hamiltonian is invariant 
under ^2 (which is isomorphic to D2), and the D3 irreducible representations reduce to those of its D2 
subgroup as Ai , E\ — )■ A (symmetric) and A2 , -Bp —> B (antisymmetric) . 

All physical operators, such as the Hamiltonian and transition operators, are expressed in terms of 
the building blocks of Eq. (H). With the help of the transposition P(12) and the cyclic permutation 
P(123) we construct a set of generators of the algebra of U{7) that transform as tensor operators under 
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the rotation group 5*0(3) and the point group chain D3 D D2 

Dp..^ ^ {p^s-s^Pp)'^^^ {Ep,B), Dx.rn ^ {p\s- s^pxi^^ (i^A,^), 

Ap,m = i{pls + s'^pp)'^^} {Ep,B), Ia,™ = i{p\s + s^px)i\'' {Ex, A), 

G% ^ {p\px+p{ppf^ {Ep,B), G^x,m ^ (pIpp ^ PiP^^)'^ {Ex, A), (5) 

GZm = {Plpp+Pipx)^n} {AuA), G^)„ = ^{plpx-pipp)ll^ (^2,5), 

fis = sh {Ai,A) , 

with I = 0, 1,2. Here Pp,„i = (— l)^^™Pp,-T?i and p^.m = (^1)^ ™PA,-m- The U{7) Hamihonian can be 
expressed in terms of scalar products of the operators of Eq. The eigenvalues and corresponding 
eigenvectors can be obtained exactly by diagonalization in an appropriate basis. 

3 Geometry and intrinsic-collective structure 

A more intuitive geometric interpretation of algebraic Hamiltonians can be obtained by using mean-field 
techniques. The main ingredient is the introduction of coherent (or intrinsic) states as variational wave 
functions which for a system of bosons have the form of a condensate of N deformed bosons Q . In the 
present case of U{7), the condensate can be parametrized as 

\N;r,x,0) = -^{bir\0), (6) 

with 



5'^ + rcosxpl + rsinx (cos^pj, .J, + sin0pj,_y) /\/l + r^ . (7) 

The geometry is chosen such that the xy plane is defined by p and A, with the a;-axis along A and the z-axis 
perpendicular to this plane. The two vectors p and A are parametrized in terms of the three Euler angles 
which are associated with the orientation of the three-body system, and three internal coordinates which 
are taken as the two lengths of the vectors rp and r^, and their relative angle 6: rx' ^p = rxrp cos . The 
two lengths are parametrized in terms of the hyperspherical radius r and the hyperangle x- fx = ?'C0Sx 
and Vp = r sin x • The hyperradius r is a measure of the size of the system, whereas the hyperangle x 
and the angle 9 determine its shape ||l0| . Since the Hamiltonian is rotationally invariant, there is no need 
to introduce the Euler angles in Eq. (0) . The expectation value of the Hamiltonian in the coherent state 
defines its energy surface 

EN{r,x,0) = {N;r,x,0\H\N;r,x,e) . (8) 

The minimum of the energy surface determines the equilibrium values (r, x,9) of the geometric variables. 
The corresponding condensate wave function \N;r,x,d) represents the equilibrium shape of the three- 
body system. 
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In order to analyze the vibrational and rotational excitations, it is convenient to split the Hamiltonian 
into an intrinsic (vibrational) and a collective (rotational) part ||ll|] 

H = H^^^ + Hc. (9) 

The intrinsic part, by definition, annihilates the equilibrium condensate and has, up to an overall constant, 
the same energy surface as the original Hamiltonian 

{N;r,x,9\H-^^^\N;r,x,0) = EN{r,x,0) - Eq . (10) 

For r > the rotational symmetry is spontaneously broken. Although in this case the condensate of 
Eq. (|^) is deformed, it is still an eigenstate of the rotationally invariant intrinsic Hamiltonian. Since H 
and i^j^^ have the same energy surface, one can extract the normal modes of the system by carrying out 
a normal mode analysis on . The collective part is the remainder of the Hamiltonian He = H — H-^^^ . 

Rather than starting from the most general algebraic Hamiltonian, analyzing its energy surface, de- 
termining its equilibrium shapes, and decomposing the Hamiltonian for each one of them into an intrinsic 
and collective part [^ , here we start from a given equilibrium shape and construct the corresponding 
intrinsic and collective Hamiltonians to analyze the vibrational and rotational excitations. In this paper 
we discuss the nonlinear equilibrium shapes of three objects, see Fig. 0. 

3.1 Intrinsic Hamiltonian 

Nonlinear configurations are characterized by the equilibrium values of the geometric variables 

f = R > 0, x = /3, ^ = 7^0,^. (11) 

The corresponding intrinsic Hamiltonian can be expressed in the form 

^int = ApIPi+BP^P2 + CpIP3 

+D (Ptp2 + p^p^) + E [pIPs + PjPi) + F {P^P^ + pIp^) , (12) 

with 

Pi = pI-pI+pI-p{^R^sU^ , 
P^ = sin^ f3pl-p{- cos^ PpI'pI , 

Pi = sin(2/3)pt .pt _cos7 (sin2/3pt .pt ^cos2/3pt .pt j . (13) 
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3.2 Collective Hamiltonian 



By construction, the collective Hamiltonian has a completely flat (or structureless) energy surface, which 
does not depend on the geometric variables (r, x,^)- Its most general form can be found by examining 
the structure of the condensate boson which defines the energy surface. The condensate boson of Eq. (^) 
is related to the scalar boson by a unitary transformation 



(14) 



with cos a = l/vT+r^ and sin a = r/vT+r^ . The operators Lpjn and F2 can be expressed in terms 
of the generators of Eq. (|) as Lp^„i = {GAlm + G^^m)/^ and F2 = -V^GaI ■ The generators of the 



transformation of Eq. (14) all belong to the 50(7) subgroup of U{7), which is generated by 



A 



J A(l) 



^A2fl ' ^ 



A2.m " 



(15) 



For this reason, the collective Hamiltonian can be expressed in terms of the two-body part of scalar 
products of these generators and has the following decomposition 



He. = H. 



H, 



(Ex 



H, 



(16) 



where 



[Ai) 



H, 



{Ep 



(Ex) _ 



Ki : Ap- Ap + Ax - Ax : +2k2 : 



Ap-Ap-Ax-Ax ■.+2k'^ 



A^-Ax+Ax- A, 



G(i).g(i)+g(i).g(i) 



(17) 



3 » ^ . ^ ^p < ^p 

Here Lm^ — v^gJ^^^ ^ represents the angular momentum operator, and :: denotes normal ordering. Note 

(2) (2) I I 

that the Ai term G\^ ■ G\^ does not appear in Eq. (17) since it is not independent 

(18) 



^(2) A(2) _ ^(1) ^(1) 1^(1) ^(1) ^(1) Ml),r)MO) MO) 
^A2 ■ '~'A2 — ' ^P ' ^\ ~ ^Ai ■ ^Ai + ^<-^A2 ■ ^A2 ■ 



The collective Hamiltonian shifts, splits and generally mixes the bands generated by ffjnt- 



4 Asymmetric tops 



We first discuss nonlinear XYZ and XY2 configurations which correspond to asymmetric tops. 
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4.1 Nonlinear XYZ configurations 

For a system of three different objects XYZ, there are no restrictions arising from the Cg point group 



symmetry The equihbrium shape is given by Eq. (11), and hence all terms in the Hamiltonian of 
Eq. are allowed. The intrinsic part is given by Eqs. (|l^) and (|l3|), and the collective part by Eqs. ( p^ ) 
and (^7|). The normal modes can be obtained by carrying out a normal mode analysis on H-y^^^. This is 
done by expressing it in terms of a deformed boson basis, which is spanned by the condensate boson of 
Eq. (I) with r = i?, X = /3 and 6* 7 (7^ 0, tt) 



bl 



/ + i?cos/3p!,_^ + i?sin/3(cos7p),_^+sin7p); ,^) /\/l + B? , (19) 



and six additional orthonormal fluctuation bosons 



hi 


= -Rs^ + 


cos/3pi_;j. +sin/3(cos7pj,^^ + sin7pj,,^) 




= sin/3p^_^ 


- cos/3 (cos7pj,^^ + sin-fpl^y) , 




= sin/?p^,^ 


+ cos (3 {sin J pI.^- cos J ply) , 


b\ 


= ' 




b\ 


= Pi,z ' 




bl 


= cos/3p1_^ 


-sin/3(sin7p), ,, -cos7p]; y) . 



(20) 

The mean-field Hamiltonian ^ H^^yit obtained by the usual Bogoliubov treatment of bosonic systems 
which amounts to replacing the condensate boson operators 6j and be by ^/N and keeping only terms to 
leading order in N |l^ 

"^-^int = eublbu + e-vblbv + ewbl^byj + euvib^by + blbu) 

+euw{bibnj + bibu) + e„„(6^5u, + bl,by) , (21) 



with 



e„ = iANR^ , 

e„ = BNR'^sm'^{2f3)/{l + R^) , 

e„ = CiVi?2 sin2(2/3)sin2 7/(1 + , 



Cuv = 2DNR^sin{2l3)/y/l + R^ , 



euw = 2ENR^ sin(2/3) sin7/ ^1 + , 

eyro = FA^i?^ sin^ (2/3) sin 7/(1 + i?2) . (22) 

The spontaneously broken rotational symmetry (i? > 0) ensures that only the vibrational bosons 5| 
(z — u, V, w) show up in ^ H:^^-^. The intrinsic modes of the system which correspond to small oscillations 
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about the miiiiniuni of the energy surfaee can be obtained by diagonahzing ^-ffj^t- The bosons 6| 
(i = 1, 2, 3) are associated with rotations of the equilibrium condensate, and hence can be identified with 
the Goldstone modes of the spontaneously broken rotational symmetry 

The mean-field or Bogoliubov image of the collective Hamiltonian has the form 

''He = Vu{bi-buf + Vv{bl-Kf+r]^{bl-b^f 

+m {bl -bif + m {b\ - b2f + m (4 - b^f 

+ 'nUV {bi - by) {bl - by) + rjyyj {b\ - by) [b]^ - by,) + r]yy, {bl " by) " by,) 

+ V12 {b\ - bl) (bl - 62) + mV {bl - 63) {bl - by) + {bl - 63) {bl - by,) , (23) 



with 



r]u = —N {ki — k'i cos(2/3) -|- k" cos7sin(2/3)} , 

T]y = -N {ki + [K2sin^ J + K4Cos'^ j]} /{I + R^) 

-N {k[ cos(2/3) + sin^ 7} /(I + R^) 

+JVKiCOS7sin(2/3)/(l + ij2) , 
rjy, = —N {ki + [k2 (cos^ 7 cos^ K4 sin^ 7 cos^(2/3)] } /(I + -R^) 

-N {k[ cos(2/3) + R'^K'2 [- sin2(2/J) + cos^ 7 cos^(2/3)] } /(I + i?^) 

+N {k'I - R^k'I cos(2/3)} cos7sin(2/3)/(l + R^) , (24) 



and 



m = - N {ki + R^ [k2 + K3 sin'^ f3]} /{I + R^) 

-N {k[ + i?2 [4 cos(2/3) -t- 24 sin^ /?] } /(I + R^) 
-NR^ [4' + 24'] cos7sin/3cos/3/(l + R^) , 

Tj2 = -7V{Ki+i?2[K2+K3COs2/3]}/(l + i?2) 

-N {-k[ + R^ [-4 cos(2/3) - 2k^ cos^ /?] } /(I + R^) 
+NR^ [k2 - 24'] cos 7 sin /3 cos /3/(l + i?^) , 
7,3 = _Ar{Ki+E2 [K2(cos27sin2(2/3)+cos2(2/?)) +K3 + K4sin27sin2(2/J)]}/(l + i?2) 

~N {-k'^cos{2I3) + R^ [4 (- cos^ (2/3) 008^7 sin^ (2/3)) - 2^ eos(2/3)] } /(H- i?^) 

-N {K'{ + R^[-Klcos{2P) + 2Kl])cos-ism{2(3)l{l + R^) . (25) 

The only nonvanishing coupling terms are r^uwi Vvw, V12, VSv and r]sy. 



r]yy = 2N{k[ sin(2/3) + k'I cos7COs(2/3)} ^ 



r]yy, = -2NK'{sm-f/\/lTR^ 
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ry^„ = -7Vi?2[K2-K4]sin(27)cos(2/3)/(l + i?2) 
-NR^k'^ sin(27) cos(2/3)/(l + R^) 

-Ari?24'sin7sin(2/3)/(l + , (26) 

and 

7yi2 = -Afi?^K3COS7sin(2/3)/(l + i?2) 
-2iVi?2^'2 cos7sin(2/3)/(l + 

~N {2k'/ - i?2 [4' cos(2/3) + 24]} /(I + i?2) , 
r;3„ = -7Vi?2[_«;2 + ^4]sin(27)sin(2/?)/(l + i?2) 

+NR^k'^ sin(27) sin(2/3)/(l + i?^) 

+N {2k'/ + i?2 [-k'2'cos(2/3) + 2k(j']} sin 7/(1 + i?^) , 
773„ = -iVi?2 [k2 - K4] sin(4/3) sin^ 7/(1 + R^) 

+N {2k'i + [24 cos(2/3)(1 + cos^ 7) + 2K^j] } sin(2/3)/(l + i?^) 

+iV{2K/cos(2/3) [_K^'cos(4/3) + 2K3Cos(2/3)]}cos7/(l + i?2) . (27) 

4.2 Nonlinear XY2 configurations 

For nonlinear XY2 configurations the relevant point group is C21, [|l2|. If we label the identical objects 
by 1 and 2, the Jacobi coordinates of Eq. (^) with mi — m2 become perpendicular. The corresponding 
equilibrium shape is 

r = R > , X = /?, ^ = 7 = 7r/2. (28) 

The permutation symmetry associated with the interchange of the two identical objects S2 (isomorphic 
to D2) is determined by the transposition P{12). This excludes from the Hamiltonian all terms that 
change the number of pp bosons by an odd number. 

Accordingly, the intrinsic Hamiltonian for nonlinear XY2 configurations is given by Eqs. (^|) and ( |l3| ) 
with the shape variables of Eq. ( ^8|) , and 

E ^ F ^ . (29) 

The normal modes are obtained as before by carrying out a mean-field analysis of the intrinsic Hamiltonian 

'^^int = ^ublbu + Cyblb^+e^bl^b^ + euviblK+blbu) . (30) 

In this case, the two radial modes, u and v, are decoupled from the angular mode w. The normal radial 
modes are obtained by diagonalizing a 2 x 2 matrix, and correspond to the symmetric and antisymmetric 
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stretching modes. This diagonaUzation is the algebraic analogue of the orthogonal transformation of 
Delves coordinates {i.e mass-scaled Jacobi coordinates) to the normal coordinates in the small displace- 
ment limit |]l3| . 

Similarly, the collective Hamiltonian is given by Eqs. (||) and (|l|) with 

H^c^"^ = . (31) 

Its Bogoliubov image has the form 

^i/c = Vu {bi - buf + Tj, {bl - + r]^ {bl - b^f 
+m {b\ - bi f + m {bl - b2f + 773 {bl - 63)' 

+yiuv {bi - bu) {bl - 6„) + 7/3«, {bl - foa) {bi - b^) . (32) 

In addition to the 7]i terms, the only coupling terms present are ?7„„ and rj^w In comparison with the 
nonlinear XYZ configuration, the coupling terms Tyuto, 'Hvw, ?7i2 and 7731, are missing. This is a consequence 
of the different point group symmetries. 



5 Oblate symmetric top 

Nonlinear X3 configurations have the equilibrium shape of an equilateral triangle, which is characterized 
by two orthogonal Jacobi coordinates of equal length. The geometric configuration has one 3-fold sym- 
metry axis, three 2-fold symmetry axes perpendicular to it and one reflection plane. The corresponding 
point group is D:^^ whose classification is equivalent to that of its subgroup 1)3 (which is isomorphic 
to the permutation group of three identical objects ^3) and parity. In comparison with nonlinear XY2 
configurations, for three identical objects the Jacobi coordinates not only are perpendicular, but also 
have equal lengths. The equilibrium shape is then characterized by 

f = i?>0, )^ = /3 = 7r/4, e ^ ^ ^ tt/2 . (33) 

The permutation symmetry associated with the interchange of the three identical objects 6*3 (isomorphic 
to 1)3) is determined by the transposition -P(12) and the cyclic permutation P(123). This excludes from 
the Hamiltonian all terms that are not scalar with respect to S'3 . The S'3 invariant Hamiltonian (or mass 
operator) is discussed in more detail in ^ . Here we briefly present the results that are relevant to the 
present article. 

The intrinsic Hamiltonian for nonlinear X3 configurations is given by Eqs. (|l|) and (|l|) with the 
shape variables of Eq. (^3|), and 

B = C , D = E = F = . (34) 
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As before, the normal modes are obtained by a mean-field analysis 

^i/ij^t = eubih, + e,{blK + blb^) . (35) 

fn this case, the two radial modes, u and v, and the angular mode w are decoupled, and correspond 
to the fundamental vibrations (see Fig. ^) of the configuration of Fig. |]. The first term represents the 
symmetric stretching mode u, whereas the second term corresponds to a degenerate doublet (e^ — £w) of 
the antisymmetric stretching mode v and the bending mode w. This is in agreement with the point-group 
classification of the fundamental vibrations of a symmetric-top X3 configuration | [T^ . 
The collective Hamiltonian is given by Eqs. (|l|) and (0) with 

H^c^"-^ ^ Hf'^ ^ , (36) 

and its Bogoliubov image has the form 

^ffc = vu{bi^buf+^,[{bi-b,f + {bi~b^f] 

+m [ib\ - b,f + [bl - b^f] + m {bl ~ b^f . (37) 

The terms which depend on the vibrational bosons b\ {i = u,v,w) cause a shift in the vibrational 
frequencies, whereas the terms involving the rotational bosons 6| {i — 1, 2, 3) are responsible for in-band 
rotational splitting. The equality of the moments of inertia about the x and y axis (ryi ~ r/2), and the 
equality of the coefficients rjy = rj^ reflect the axial symmetry of the underlying shape. In comparison 
with the nonlinear XYZ and XY2 configurations all coupling terms rjij are missing. 

The above analysis presents a clear interpretation of the various interaction terms in the intrinsic and 
collective Hamiltonians. The two terms in i^j^^ proportional to A and B ^ C contribute only to the 
eigenfrequencies of the vibrational modes. The K3 and K4 terms in He contribute only to the rotational 
modes, whereas the ki and K2 terms both effect the vibrational and the rotational modes, and hence 
represent vibration-rotation couplings. The equality of the frequencies of the antisymmetric stretching 
and the bending modes, and the equality of the moments of inertia about the x and y axis, show that the 
above Hamiltonians correspond to an oblate symmetric top with the threefold symmetry axis along the 
z-axis. In Q it was shown that U{7) provides a complete classification of the vibrational and rotational 
excitations of the oblate symmetric top. 

6 Summary and conclusions 

We have presented an algebraic treatment of three-body problems. The relative motion of the three- 
body system is treated by the method of bosonic quantization, which for the two relative vectors gives 
rise to a U{7) spectrum generating algebra. The model space is spanned by the symmetric irreducible 
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representation [N] of U{7). In particular, we have considered nonlinear XYZ, XY2 and X3 configurations. 
The U{7) interaction terms have been interpreted geometrically by means of coherent states and a normal 
mode analysis. We have studied the effect of the permutation symmetry among the identical particles on 
the coupling between the vibrational and rotational modes. For XYZ and XY2 the geometry is that of 
an asymmetric top, whereas for that of an oblate top. The ensuing algebraic treatment of the oblate 
top has found useful applications both in molecular physics (X3 molecules |6j) and hadronic physics 
(nonstrange qqq baryons 

In conclusion, we have shown that U{7) which was previously introduced to describe X3 configurations 
^, also provides a spectrum generating algebra for XY2 and XYZ configurations. The latter are 
important to describe nonlinear molecules and strange baryons. In this respect, U (7) provides a unified 
treatment of vibrational and rotational excitations of both symmetric and asymmetric tops. 
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Figure 1: Geometry of a three-body system. 




P 

u vibration v vibration w vibration 



Figure 2: Schematic representation of the normal vibrations of a nonlinear X3 configuration. The Jacobi coor- 
dinates are indicated by the dotted lines. 
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